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ABSTRACT 

We build the Yunnan-III evolutionary population synthesis (EPS) models by using 
the MESA stellar evolution code, BaSeL stellar spectra library and the initial mass 
functions (IMFs) of Kroupa and Salpcter, and present colours and integrated spectral 
energy distributions (ISEDs) of solar-metallicity stellar populations (SPs) in the range 
of lMyr-15 Gyr. The main characteristic of the Yunnan-III EPS models is the usage of 
a set of self-consistent solar-metallicity stellar evolutionary tracks (the masses of stars 
are from 0.1 to lOOM©). This set of tracks is obtained by using the state-of-the-art 
MESA code. 

MESA code can evolve stellar models through thermally pulsing asymptotic giant 
branch (TP-AGB) phase for low- and intermediate-mass stars. By comparisons, we 
confirm that the inclusion of TP-AGB stars make the V — K, V — J and V — R colours 
of SPs redder and the infrared flux larger at ages log(t/yr) ^7.6 (the differences reach 
the maximum at log(i/yr)~8.6, ~0. 5-0. 2 mag for colours, ^ 2 times for A'— band flux). 
We also find that the colour-evolution trends of Model with-TPAGB at intermediate 
and large ages are similar to those from the starbvrst99 code, which employs the 
Padova-AGB stellar library, BaSeL spectral library and the Kroupa IMF. At last, we 
compare the colours with the other EPS models comprising TP-AGB stars (such as 
CB07, M05, VIO and POPSTAR), and find that the B - F colour agrees with each 
other but the V — K colour exists larger discrepancy among these EPS models (~1 mag 
when 8<log(t/yr) <9). 

The stellar evolutionary tracks, isochrones, colours and ISEDs can 
be obtained on request from the first author or from our website 
dhttp : / / wwwl .ynao.ac.cn/ '^zhangfh /) . Using the isochrones, you can build your 
EPS models. Now the format of stellar evolutionary tracks is the same as that in 
the STARBi]RST99 code, you can put them into the starburst99 code and get the 
SP's results. Moreover, the colours involving other passbands or on other systems 
(for example, HST F439W^ — colour on AB system) can also be obtained on 

request. 

Key words: stars: evolution — stars: AGB and post-AGB — galaxies: stellar content 
— galaxies: fundamental parameters — galaxies: evolution 



1 INTRODUCTION 

Evolutionary population syn thesis (EPS) method was first 
introduced bv iTinslevI ()l968h . and mainly is used to study 
the unsolved stellar population (SP). Key ingredients in EPS 
models are: (i) the library of evolutionary tracks used to 
calculate isochrones in the colour-magnitude diagram; (ii) 
the library of stellar spectra adopted to derive the inte- 
grated spectral energy distribution (ISED) or magnitudes 
and colours in suitable passbands; and (iii) the initial mass 
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function (IMF) used to evaluate the relative proportions of 
stars in the various evolutionary phases. 

Stellar evolutionary tracks and isochrones play an im- 
portant role in EPS models. However, many EPS mod- 
els use the combination of seve ral stellar evolu tionary li- 
braries [such as the models of IWorthevI (|l994l '). PEGASE 
(IFioc fc Rocca- Volmerangel 1 19971 . 1 19991 ) ." etc.]. This implies 
that the basic stellar properties are dis continuous along 
the composite tracks (jMarigo et al.l [20081 . hereafter M08). 
Moreover, due to theoretical and observational limitations, 
many stellar spectra libraries only cover the limited param- 
eter (such as effective temperature T^s, surface gravity logg. 
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metallicity [Fe/H], etc.) ranges. Therefore, some EPS models 
use the combination of several stellar spectra libraries. 

Based on the above-mentioned facts, in this paper we 
will build the Yunnan-III EPS models and present the 
ISEDs and colours of solar- metallicity SPs by using a set 
of self-consistent stellar evolution models (computed from 
Modules for Experiments in Stellar Astrophysics [mbsa] 
code, its properties will be described in Section [ 2]), a stellar 
spectra library (BaSeL-2.2, iLeieune et al.lll997l.1l998l . here- 
after LCB97, LCB9 8) and th e IMFs of lKroupa et'aiTl|200ll . 
hereafter KOI) and ISalpeteJ (|l955l . hereafter S55). Before 
the construction of the Yunnan-III models, we have built 
two kinds of EPS models. The Yunnan-I mo dels are con- 
struc ted for SPs without binary interactions (jZhang et al.l 
|2002|) by using the rapid single stellar evolution algorithms 
(|Hurlev et al.ll2000l . hereafter HPTOO). The Yunnan-II mod- 

eijr 

are cons tructed for both SPs wi th and without binary 
interactions (jZhang et al.ll2004l . |2005| ) by using Monte Carlo 
simulation and the r apid single and bin ary stellar evolution 
algorithms CHPTOO. iHurlev et al.|[2002l . hereafter HTP02). 

Another main reason of building the Yunnan-III models 
is that stars on thermally pulsing asymptotic giant branch 
(TP-AGB) pha se are often neg lected or added following indi- 
vidual recipes (|Marastonll2005l . hereafter M05) in EPS mod- 
els. This is because computations of the TP-AGB phase by 
means of complete evolution code are very demand ing in 
terms of computational time (|Marigo fc Girardill2007l . here- 
after MG07). So far, there exists several EPS models that 
have considered the contribution of TP-AGB stars. These 
models would be recalled in Section|3l According to the used 
technique of dealing with TP-AGB stars, these EPS models 
are divided into two groups. 

With 'fuel consumption' technique, iMarastonI (|l998l . 
hereafter M98) and M05 have considered the contribution 
of TP-AGB stars in her EPS models. The advantages of 
the 'fuel consumption' approach are summarized by M05. 
Based on her models, she and her colleagues have investi- 
gat ed the redshift z ~ 2 pass ive and star- forming galax- 
ies (iMaraston et al.ll2006l . I2OI0I ) and luminous red galaxies 
l|Maraston et al.ll2009l ) 

With 'isochrone synthesis' technique, some models 
[such as the mode ls of CBll and their older CB07 ver- 
sion (|Bruzuail I2OII ), Flexible Stella r Population Synthesis 
(hereafter FSPS. ICoiiro v ct al."2009': 'Coii rov fc Gunnll2010l : 
IConrov et aLll2010l ), IVazdckis ct al (201ij, hereafter VIO) 
STARBURST99 (hereafter S B99, iLeitherer et al.l 19991, 2010l: 
Vazquez fc Leithereill2005l)-v6.02, POPSTAR jMoUa etlo] 



20091 ). Yunnan-I and II l|Zhang et aLll2002l . [iooa ). etc.] have 
included TP-AGB stars. In these models, the evolution of 
TP-AGB stars is obtained either by perturbation function 
or by a synthetic code. 

The Yunnan-I and II EPS models have used TP-AGB 
evolution from perturbation function of the HPTOO rapid 
stellar evolution code. HPTOO claimed that this method can 
produce the core-mass (Mc) and luminosity (L) relation of 
the third dredge- up (TDU). 

Then, what is a TP-AGB synthetic code? In a TP- 
AGB synthetic code, the stellar evolution is described by 
means of simplified relations (such as Mc-L relation). It 
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can be divided into detailed and simple codes, (i) In a de- 
tailedly synthetic code, the simplified relations are derived 
from complete stellar models, and convective dredge-up and 
mass loss are tuned by means of a few adjustable param- 
eters (M G07). A very detailedly s ynthe tic code was d evel- 
oped by iMarigo et all (| 19961 . Il998l ) and iMarigol (|200ll , and 
references therein) . MG07 and their colleagues have used an 
explicitly synthetic code, in which the TDU and hot bot- 
tom burning (HBB) are considered, to present the evolution 
of TP- AGB stars and the firs t isochrones considering the 
TDU iMarigo fc Girardil [iooH ) . (u) In a simply synthetic 
code, the TDU and the over-luminosit y above the M c -L re- 
lation caused by HBB are neglected. Bertelli et al.l (|l994l . 
hereafter BBCFN-1994) and iGirardi et al.l (|200G| . hereafter 
GOO) have used a simply synthetic method to present the 
TP-AGB evolution (M08). 

Among those 'isochrone synthesis' EPS models includ- 
ing TP-AGB stars by a synthetic code, the TP-AGB evolu- 
tion of the GBll, GB07 and FSPS models is from an explic- 
itly synthetic code, while that of the VIO and POPSTAR 
models is from a very simply synthetic code. Looking at 
the ingredients of the above EPS models including TP-AGB 
stars by a synthetic code (refer to Section 3), we find that 
all the above-mentioned models used the TP-AGB evolu- 
tionary tracks from Padova group. The reason of the above 
mentioned situation is that the use of the so-called synthetic 
code is the only viable alternative at that time in order to 
provide extended grids of TP-AGB models that reproduce 
basic observational constraints (MG07). 

However, it is no doubt that stellar evolution model 
is another method of obtaining TP-AGB evolution. Stellar 
evolution models often parameterize physical proc esses, this 
is contrary to the synthetic method (|Herwig|[2005l ). 

In this work, building the Yunnan-III models with 
'isochrone synthesis' technique, we will compute the stel- 
lar evolutionary tracks including TP-AGB phase by the 
MESA code, which can compute the evolution of low- 
and intermediate-mass stars from zero-age main sequence 
through TP-AGB phase. Therefore, the Yunnan-III EPS 
models can be compared with those EPS models that have 
considered TP-AGB stars, and can be an alternative one in 
galaxy studies. It is worth noting that the TDU is consid- 
ered in the MESA star and massive TP-AGB stars show the 
expected HBB behavior, including the avoidance of the C- 
star phase for a 5Mb, ^=0.01 st ellar model track despite 
eflicient TDU (|Paxton et al.ll201ll ). 

The outline of the paper is as follows. In Section 2 we 
describe the Yunnan-III EPS models. In Section 3 we re- 
call some EPS models considering TP-AGB stars. In Sec- 
tions 4 and 5 we present the results (evolutionary tracks, 
isochrones, colours and ISEDs of SPs) and comparisons. Fi- 
nally we present a summary and conclusions in Section 6. 



2 MODEL DESCRIPTION 

In this section, we describe the employed ingredients (includ- 
ing the MESA stellar evolution code, BaSeL stellar spectra 
library and IMFs) in the Yunnan-III EPS models. The main 
characteristics are listed in Table [1] 
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2.1 MESA stellar evolution code: input physics and 
parameters 

Stellar evolution models form the basic of EPS studies, and 
therefore we devote much of the remainder of this section to 
describing aspects of the MESA models in details. 

The MESA s tellar evolution code was built by 
IPaxton et alj (|201lh , and beg an as an effort to improve upo n 
the EZ stellar evolution code (|Eggletonlll97ll : |Paxtonll2004h . 
During the construction, its design and implementation were 
influenced by a number of stellar evolution and hydrody- 
namic codes (therefore, it can calculate the evolutions of 
massive stars). MESA is a state-of-the-art stellar evolution 
code, it uses the most updated equation of state (EOS, can 
handle the low mass and high degree of degeneracy), opac- 
ity, nuclear reaction rates, element diffusion data, and atmo- 
sphere boundary conditions. Moreover, MESA code employs 
modern software engineering tools and techniques to target 
modern computer architectures, and its numerical and com- 
putational methods employed allow it to consistently evolve 
stellar models through challenging phase, for example, MESA 
code allows non-stop computations through the helium core 
flash toward the end of the AGB evolution for low-mass 
stars (jPaxton et al At last, MESA code uses the adap- 

tive mesh refinement and sophisticated time-step controls, 
therefore, it is able to solve the fully coupled structure and 
composition equations simultaneously. In the following, we 
describe the input physics and parameters used by us in the 
MESA code. 

Network: the 'basic' network of eight isotopes (^H, 
^He, "^He, ^^C, "N, 1*^0, ^°Ne and ^^Mg) is the defauh value. 
We use the 21 isotope reaction netwo rk (approx21.net), in - 
spired by the 19 isotopy networlfl in IWeaver et all (|l978h . 
that is capable of efficiently generating accurate nuclear en- 
ergy generation rates from hydrogen burning through silicon 
burning. This network includes linkages for PP-I, steady- 
state CNO cycles, a standard a— chain, heavy ion reactions, 
a nd aspects of photo disintegration in ^'*Fe (see Section 7.3 
of IPaxton et allboill ') . 

EOS: we u se the MESA, HELM jTiinmes fc Swestvl 

I2OOOI ) and PC (|Potekhin fc ChabrieJ I2OI0I ) EOS tables 
The MESA EOS ta ble is constructed from t he 2005 up- 
date of the OPAL jRogers fc Navfonovl |2002| ) and SCVH 
(jSaumon et al.ll 19951 ) EOS tables. 

Opacity: we use the opaci ty of Grevesse fc Sauvaj 

1 19981 . hereafter GS98) instead of iGrevesse fc Noeld jimS . 
hereafter GN93, by default), with t he low temperature opac- 
ities taken from iFerguson et al.l (|2005|) with metal ratios 
given by GS98. The low T opacities of lFerguson et all (|2005l ') 
include the effects of molecules and grains on the radiative 
opacity. 

Mixing length parameter: for M < SOM©, mixing 
length parameter a = 2; for M > SOM©, a = 1. 

Overshooting: MESA treats convective mixing as a 
diffusive process with a diffusion coefficient. The overshoot 
mixing diffusion coefficient Dov = -Dconv,oexp( — j^^), in 
which -Dconv.o is the mixing-length-theory derived diffusion 
coefficient at a user-defined location near the Schwarzschild 



boundary, Ap,o is the pressure scale height at that location, 
z is the distance in the radiative layer away from that lo- 
cation, and / is an adjustable parameter. By default, the 
overshooting is turn off. We switch on it when star mass 
M > I.8M0, and turn off it when M < I.IMq. In the 
mass range of 1.1 — 1.8 Mq, overshooting is gradually en- 
abled. Moreover, convective overshooting is only allowed at 
locations with m > 0.001 Af (by default). At the upper and 
lower convective boundaries for H-burning and non-burning 
convective zones, the adjust parameter / = 0.01. 

Mass loss: two kinds of mass loss are needed to be 
mentioned. The first one is the mass loss of red giant 
branch (RGB) and AGB stars. B y defaul t, ther e is no 
RGB and AGB wind. We use the iReimersI (Il975l, Mr = 



4 X 10""77rL7?/M,77r = 0.5) and |BloeckerJ (| 19951 . Mb 



2 similar to Frank Timmes' APPROX19, sec MESA/data/net_ 
data/nets/README 



4.83 X IQ-^tibL'^-^ M-^-^Mn/m, m = 0.1) mass loss laws 
for RGB and AGB stars, respectively. The second one is the 
enhanced mass loss due to rotation. By default, it is consid- 
ered and the coefficient is 0.43. In this work, we set it for 
the reason of comparison with the other EPS models. 

Atmosphere: the default atmosphere is a simple pho- 
tosphere model, i.e. do not integrate and just estimate for 
T = 2/3. We use the default one. 

Comments: (i) the TDU: in the MESA code, the 
TDU is considered and massive TP-AGB stars show the 
HBB behavior, (ii) the set of input parameters: 
using the above set of input parameters and physics, 
we obtain the luminosity of 1 M© star at an age of 
4.57x10^1-, log(L/L0)=-O.Ol for the GS98 opacity while 
log(Z//LQ)=0.008 for the GN93 opacity. These luminosities 
are within one hundredth of a magnitude. 



2.2 BaSeL stellar spectra library 

The BaSeL-2.2 stellar spectra library (LCB97, LCB98) 
is a semi-empirical one and provides an extensive and 
homogeneous grid of low-resolution theoretical flux dis- 
tributions in the range of 9.1— 160 000 nm and synthetic 
UBVRIJHKLM colours for a large range of stellar pa- 
rameters: 2000 < Toff/K < 50 000, -1.02 < log 5 < 5.50 and 
-5.0 < [Fe/H] < 1.0. We adopt the BaSeL stellar spectra 
library for the following reasons, (i) The BaSeL library has 
wider cover age of wavele ngth in comparison with other li- 
braries (see lBruzualll201 ih . so the EPS models based on this 
library can be used to the studies of galaxies at different 
redshift via multi-band SED fitting method, (ii) The BaSeL 
library provides wider coverage of metallicity in comparison 
with the empirical stellar spectra libraries. Our next work 
is to give the EPS models at the other metallicities (even at 
zero metallicity). Therefore, we use this library for the sake 
of homogeneity. 

At solar metallicity, the BaSeL library covers the region 
occupied by the TP-AGB stars on the log p-logTcff plane, 
but does not include TP-AGB stars. Therefore, we use the 
normal M-giant spectra, provided by the BaSeL library, for 
TP-AGB stars in this work. 

The situation, whether or not the adopted stellar spec- 
tra library includes the spectra of TP-AGB stars, would af- 
fect significantly the infrared (IR) results of SPs compris- 
ing TP-AGB stars. Recently, some EPS models have used 
the empirical library comprising TP-AGB stars (such as the 
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M05 model). Therefore, in the future, we will use these em- 
pirical library comprising TP-AGB stars in our models. 



2.3 IMF 

We use the IMF of KOI: 

r M-°-^°, 0.01 < M < 0.08, 



0(^)koi = < M~'^-^°, 0.08 < M < 0.50, 
I ]Vr^-''°, 0.50 < M < 100, 



(1) 



where M is the stellar mass in units of M0 . 

In order to compare with the other EPS models, we also 
use the S55 IMF: 

(p{M)s^^ ^ M-^ '^^ ,0.1 < M < 100., (2) 
where M is the stellar mass in units of M© . 

3 OVERVIEW OF EPS MODELS 

In this section, we will overview those EPS models which 
have considered the contribution of TP-AGB stars. The 
main characteristics, including stellar evolution library, stel- 
lar spectra library and IMF, are summarized in Table [T] 

To distinguish which EPS models include TP-AGB evo- 
lution by an explicitly or a simply synthetic codes, we first 
give a brief description of the Padova stellar evolutionary 
tracks and isochrones. 

(i) B BCFN-1994 ( Z=0.0 004, 0.004, 0.008, 0.02 and 
0.05) and lGirardi et al.l ij 19961 . hereafter G96, Z^IO""*) pre- 
sented the isochrones comprising TP-AGB evolution from a 
simply synthetic code. These two sets of isochrones are de- 
rived from the so-called Padova- 1994 tracks, which are used 
in the lBruzual fc CharlotI (|2003l . hereafter BC03) and POP- 
STAR models. BBCFN-1994 isochrones are the basis of the 
following two sets of isochrones. 

(ii) GOO used a simply synthetic prescription to include 
TP-AGB regime to the point of complete envelope ejection 
for stars within O.15<Af/M0 <7. In the work of GOO, the 
updated opacities and EOS, and a moderate amount of con- 
vective overshoot are used, while TDU is not taken into 
account. Also, they built the isochrones (Z=0.0004, 0.001, 
0.004, 0.008, 0.019 and 0.03, 10''-*-10^° Vr)- This set of 
isochrones is used by the VIO models. Further, combining 
with the work of BBCFN-1994, they built the Girardi-2002 
isochrones (= [BBCFN-1994] + [GOO] + [simplified TP- 
AGB]). 

(iii) MG07 presented TP-AGB evolution by using a de- 
tailedly synthetic code and some recipes (semi-empirical pre- 
scriptions). They started from the ph ysical conditio ns at 
the first thermal pulse from GOO and ICirardil ()200lh and 
computed the evolution up to the complete loss of the en- 
velope via stellar winds (M08). In their models, several 
important theoretical improvements over previous calcula- 
tions (such as, TDU and HBB) have been included. More- 
over, the prescriptions have been calibrated against near- 
IR data from the Large Magellanic Clouds (LMC) and 
Small Magellanic Clouds (SMC), including carbon star lu- 
minosity function and TP-AG B lifetime (star count) in MC 
clusters (jConrov et al.ll2009l l. M08 presented the Marigo- 
2008 isochrones based on the work of MG07 {= [BBCFN- 
1994] + [GOO] + [MG07]). Based on the Marigo-2008 



isochrones, ICirardi et al.l (|2010l . hereafter GIO) further built 
the isochrones with two kinds of corrections for low-mass 
and low-Z AGB tracks. The TP-AGB evolution of MG07 is 
used by the CBll, CB07, FSPS and SB99-v6.02 models. 

3.1 M05 and the latest version MSll 

Maraston have used the 'fuel consumption theorem' to build 
a scries o f EPS models. In th e M05 and the latest version 
(Marasto n fc Strombac^l201ll . hereafter MSll) EPS mod- 
els, TP-AGB stars are semi-empirically included and cali- 
brated with observational data. 

In the M05 models, several stellar evolutionary li- 
brariejfl, the BaSeL-3 . 1 ste llar spectra library (LCB97; 
LCB98; IWestera et all 120021) . and the S55 (0(M)s55 = 
M~°',a = 2.35) and KOI IMFs were us ed. For the spectra of 
TP-A GB stars, the empirical library bv lLancon fc Mouhcind 
(|2002l . hereafter LM02) was used. 

The MSll models were built on the basis of 
M05 models by replacing the BaSeL-3.1 library by 
four empirical h igh-to- interme diate spectra l resolution li- 
brari es, namely iPicklesI (Il998f). ELODIE (jPrugniel et all 
l2007l). STELIB jLe Borgne et all |2003| ) and MILES 
ISanchez-Blazguez et al.ll2006l ) . In addition to the above up- 
dation, they als o presented the results based on the IMF of 
IChabrieJ (|2003l . hereafter Cha03). 



In the M05 models, six metallicities {Z=10 



10- 



0.01, 0.02, 0.04 and 0.07) are included in the age range from 
3 xlO"''Gyr to 15Gyr (67 ages, 16 ages within l-15Gyr 
for ^=10"" and 0.07). 

In the MSll models, several metallicities and ages are 
included. For different stellar spectra library, the metallicity 
and age coverage is differenlQ. 



3.2 CBll and their preliminary version CB07 

Chariot fc Bruzual have built a series of EPS models. In the 
CBll and their preliminary version CB07 models, TP-AGB 
stars have been included. 

The main characteristics of the CBll models are 
as follows. T hey were buil t by u sing (i) the evolution- 
ary tracks of iBertelli et all (|2008l . hereafter BGMN08)E], 
the TP-AGB prescriptions of MG07 and GIO, (u) the 
MILES stellar spectra library, complemented at the hot 



3 Cassisi + Geneva for 10"^ < Z < 0.04 (Geneva for 10"^ < t < 
30Myr); Cassisi for Z = 10'^ and Padova for Z =0.07. In the 
above d escriptions, 'Cassisi' mea ns the work of ICassisi et al] 
lll997al) . [Cassisi et al.l lll997bl) andlCassisi et al.l l l2000h. 'Geneva' 
includ es the work of ISchaher et al. I lll992h and iMevnet et alj 
lll994l) . 

For Pickles: Z=0.02; 2.5Myr-12/15Gyr. For STELIB: Z=0.01, 
0.02 and 0.04; 200/30/400Myr-12/15Gyr. For MILES: Z = 
10-*, 10-3, 0.01, 0.02 and 0.04; 5G/2G/55M/~6M(3G)/100Myr 
-15Gyr. For ELODIE: Z = IQ-^, 10'^, 0.02 and 0.04; 6G/55M/ 
3M/100Myr - ll/15Gyr. 

^ In fact, CBll used the tracks up to 15 Mq from the models 
with updated input physics by BGMN08. For stars more mas- 
sive than 15 Mq, in the range from 20 to 120 M©, CBll used 
the so-called Padova-1994 tracks. The TP-AGB evolution of low- 
and intermediate-mass stars is followed according to the semi- 
empirical proscription of GIO. 
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Table 1. Main characteristics of EPS (M05, MSll, CBll, CB07, FSPS, VIO, SB99, POPSTAR, YUNNAN-I/II an III) models. '/' means 
that the corresponding EPS ingredient can be chosen. '+' means that the EPS ingredient is from several sources. 



Name 


stellar 


evolution library /isochrone 




stellar spectra library 


IMF 


M05 


Cassisi 


+ Geneva + Padova 




BaScL-3.1-|-LM02 


S55/K01 


MSll 








Picklea/ELODIE/STELIB/MILES + LMQ2 


S55/K01/Cha03 


CB07 


Padova 


-1994/Padova-2000 + MG07 




BaScL-3.1/STELIB/Pickles 


S55/Clia03 


CBll 


Padova 


-1994 + BGMN08 + MG07 + GIO 




MILES-l-othcr 


Cha03 


FSPS 


Baraffc 


ot al.(199S) 4- MG07 + M08 




BaScL-3.1-|-LM02 


van Dokkum (2008) 


VIO 


G96(Z 


= 10"'^) + G00(Z>10~'^) isochrones 




MILES 


unimodal/bimodal/KOl/KOl-cor 


SB99 


Padova 


-AGB/Padova-STD/Goncva-STD/Gcnc 


va-HIG 




change 


POPSTAR 


Padova 


-1994+BGS98 




BaScL-|- other 


S55/Fcrrini/K01/Cha03 


Yunnan-I/II 


Yunnar 


-I/II 




BaScL/HRES/BLUERED 


MS79/S55 


Yunnan-III 


Yunnai 


-III 




BaScL 


S55/K01 



effective te mperature end by theoretical model atmo - 



spheres fromlLanz fc Hubenvl iMartins et al 



Lanz fc Hubenv 



12005). 



( 2007f) 



iRodn'guez-Merino et al. ( 20051 ) and 

and at the cool en d by the IRTF l ibrary ( Rayner et al. l2009l l 
and the models of lAringer et al.l (|2009l ). and (iii) the Cha03 
IMF with the mass limits of 0.1 and 100 Mq. Here, it is 
emphasized that BGMN08 also presented TP-AGB tracks, 
based on the MG07 prescriptions but extrapolated to differ- 
ent c hemical compositions of the stellar envelope (jBruzuall 

The CB07 models are similar to the BC03 models, ex- 
cept the use of the MG07 prescriptions to describe TP-AGB 
evolution. Moreover, CB07 models complemented the re- 
sults at Z=Q.l. The CB07 models used (i) two sets of stellar 
evolutionary tracks (Padova- 1994 & Padova-2000), (ii) two 
forms of IMF (S55 with a = -2.35 & Cha03), and (iii) high- 
(STELIB & Pickles) and low-resolution (BaSeL-3.1) stellar 
spectra libraries. For both IMFs, the lower and upper mass 
limits are 0.1 and 100. M©. 

The CBll models have not been released, while the 
BC07 models are available. For the CB07 models, seven 
metallicities (^=10"*, 4x10'*, 4xl0~^ 8xl0'^ 0.02, 0.05 
and 0.1) and 220 ages (from lO'^ to 10^"'^ yr) are included. 



3.3 FSPS: 

FSPS models also have considered the contribution of TP- 
AGB stars. Thes e models mak e use o f (i) the non-evolving 
stellar models of iBaraffe et al.l (|l998l ) for stars in the mass 
range 0.10 < M/Mq < 0.15 and the models of MG07 and 
M08, (ii) the BaSeL-3.1 library, complemented by the av- 
erage TP-AGB spectra compiled by LM02 from more than 
100 o ptical/near-IR spectra presented in iLancon fc WoodI 
l|2000l ) (this is sim ilar to that o f M05 ), and (iii) the IMF 
form advocated by Ivan DokkumI (|2008l ): 



dn 
dluM 



(A/ > 



(3) 



in which Ai = 0.140, = 25, x = 1.3, a = 0.69 and Ah 
= 0.158. Variation of the IMF is incorporated in the char- 
acteristic mass mc. The lower and upper mass limits of 0.1 
and 100 M© are used, respectively. 



3.4 VIO: 

Vazdekis also has built some EPS models. His latest mod- 
els, VIO, make use of (i) the updated version of the mod- 
els published in G96 for Z—W~'^ (the calculations are now 



compatible with GOO) and the solar-scaled theoretical GOO 
isochrones for Z > 10"*, (ii) the MILES stellar spectra li- 
brary, and (iii) four forms of IMF: unimodal, bimodal, KOl- 
universal and KOl-revised (hereafter KOl-cor). The slope 
of IMF for unimodal and bimodal forms is in the range 
of 0.3—3.3. The lower and upper mass limits are 0.1 and 
100 M0 , respectively. 

VIO models provide the results at seven metallici- 



ties {Z= 10'", 4x10"*, 



10" 



4x10'^ 8x10"^ 0.019 and 



0.03), the age range is different at different Z: 10— 18Gyr 
at Z=10"* (only for IMF slopes < 1.8); 0.07-18 Gyr at 
Z=4xl0"*; and 0.06-18 Gyr at 10"^ <Z< 0.03. 



3.5 SB99 v6.02: 

SB99 models were built by Leitherer and his colleagues. We 
use the 6.02 version. The desc ri ption of the input physics is 
give n bvlLeitherer et all (Il999l) , IVazguez fc Leithereij (120051 ) 
and iLeitherer et^T ( 2010l ). The main characteristics are as 
follows. 



• Four sets of stellar evolutionary tracks are provided, 
with each set including five metallicities. These four sets 
are the Geneva tracks with standard and high mass- 
loss rates (hereafter Geneva-STD fc Geneva-HIG), the 
Padova tracks as updated by GOO (hereafter Padova-STD), 
and the Padova tracks with the inclusion of TP-AGB 
stars (hereafter Pad o va-AG B) following the prescription of 
IVassiliadis fc WoodI (|l993t l. For Geneva tracks, Z=10 , 
0.004, 0.008, 0.02 and 0.04. For Padova tracks, Z= 0.0004, 
0.004, 0.008, 0.02 and 0.05. 

• Five sets of stellar atmosphere libraries are provided to 
choose. BaSeL library is one of five. 

• By default, the Kroupa IMF with two mass intervals 
(the exponent a = [—1.3, —2.3], the mass boundary Mcut = 
[0.1, 0.5, 100.] Mq) is used. The IMF form can be changed 
according to the need. 



3.6 POPSTAR: 



The POPSTAR models were built by iMoUa et"aLl (|2009l) . 

and the ingredients are described as follows. 

• POPSTAR models a dopt stellar evolution ary tracks 
from lBressan et al.l (|l993l ). lFagotto et all (|l994al lbl) and G96 
(i.e. the Padova-1994 tracks) , and combine a revision of the 
Padova (iBressan et al ] Il998l , here after BGS98) isochrones 
used in iGarcia- Vargas et al.l (|l998l ). 
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• POPSTAR models mainly use the BaSeL stellar atmo- 
sphere models, which are for normal stars. For O, B and 
Wolf-Raylet (WR) stars, POPSTAR models take the non- 
loc al thermodynamic equilibrium (NLTE) blanketed models 
bv lSmithet al] (|2002l ) for Z = 0.001, 0.004, 0.008, 0.02 and 
0.04. For post- AGB and plane tary-nebula stars, POPSTAR 
models use the iRauchI (|2003l ) NLTE models, covering the 
temperature range of 50 000-190 000 K and \ogg of 5.00-8.00. 
For higher temperatures, POPSTAR models use blackbod- 
ies. These models include all elements from H to Ni, and 
they are available for two values of metallicities: Z — 0.002 
and 0.02. POPSTAR models use the first metallicity spec- 
tra for their three metal-poor {Z < 0.004) isochrones, and 
the solar abundance spectra for the three other metallicity 
{Z > 0.008) isochrones. 

• POPSTAR models include six IMFs. They are the S55 
law with a = - 2.35 in the range of 0.85-120, 0.15-100 and 
1.00-100 Mq- the lFerrini et al.l(|l990D . KOI and Cha03 IMFs 
with masses between 0.15 and 100 M0. 

At last, POPSTAR models provide the results at six 
metallicities (10'"^, 4x10""*, 4xl0"^ 8xl0"^ 0.02 and 0.05) 
and ~106 ages (10^- - 10^°-^Vr)- 

3.7 Yunnan-I and II: 

In the Yunnan-I and II models, rapid single and binary stel- 
lar evolution codes (HPTOO, HTP02) are used. In the HPTOO 
and IITP02 codes, TP-AGB stars are obtained by perturba- 
tion function. In the Yunnan-I models, only for SPs without 
binary interactions, t he BaSeL-2.2 and se veral IMFs [S55 
with different a and iKroupa et al] (| 19931 . hereafter K93)] 
are used. In the Yunnan-II models, for both SPs with and 
without binary interactions, several [including the BaSeL- 
2.2, BLUERED jBcrtcU i ct al , .200a . ~ 0.3 A), and HRES 
(jGonzalez Delgado et al. I booiT ~ 0.1 A)] stellar spectra li- 
braries are used. For SPs with binary interactions, both 
IMFs for the primaries and secondari es are needed. We use 
the IMFs of S55 with a = -2.35 and lMiller fc Scald (|l979l ) 
for the primaries. The lower and upper mass limits are 0.1 
and 100 Mq , respectively. 

Both models can provide the results at seven metallic- 
ities {Z=W~*, 3x10"'', 10-^ 4xl0-^ 0.01, 0.02 and 0.03) 
and in the range from 10^ to lO^'^^yr. 



4 RESULTS AND COMPARISONS WITH SB99 
MODELS 

In this section, we will present the colours and ISEDs for 
SPs. The ages are in the range of lMyr-15Gyr. To in- 
vestigate the effect of IMF, we perform two sets of cal- 
culations: one uses the KOI IMF whereas another uses 
the S55 IMF. The magnitudes and colours for the mod- 
els with the KOI and S55 IMFs are presented in Ta- 
bles lAll and IA2I respectively. The ISEDs and isochrones 
can be obtained on request from the first author or from 
our website (http://wwwl.ynao.ac.cn/~zhangfh/). More- 
over, m Tables and IX2I we only give several Vega colours 
on Johnson-Cousin system, the colours involving other pass- 
bands (such as b — y colour) or those on AB system can also 
be obtained on request. 
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Figure 1. Stellar evolutionary tracks (grey line) and isochrones 
(red line). The masses of stars are in the range of 0.1— 100 Mq 
and at a logarithmic interval of 0.05. The ages of isochrones are 
log(t/yr) = 6, 7, 8, 9 and 10. The TP-AGB tracks are represented 
by blue symbols. 



T 
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Figure 2. TP-AGB phase on the HR diagram for stars with 1.0 
< M/Mq < 7.9 [log(M/M0) = (i-l)*0.05, i=l-20. Black points 
are for i=l-5 (I.O-I.SSMq), red are for i=6-10 (1.78-2.82Mq), 
green for j=ll-15 (3.16-5. OIMq) and blue for i=16-19 (5.62- 
7.94M0)]. 



At last, we will give the comparisons in the colours be- 
tween our and the SB99 models. The comparisons with the 
other literatures will be given in the next section. The rea- 
son, we put the comparisons with the SB99 models in this 
section, is that the SB99 code comprises two sets of stellar 
evolutionary tracks: Padova-STD and Padova-AGB (com- 
prising and neglecting TP-AGB stars). This can be used to 
compare the effect of TP-AGB stars on the results between 
ours and theirs. 
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log(r^„/K) log(7'^„/K) 

Figure 4. Comparisons between the model (green line) and interpolated (red open circles) tracks for a low-mass star [M = 1.78 M0, 
left panel) and an intermediate-mass star (A/ = 5.62 Mq, right panel). In each panel, the tracks of the two adjacent stars (black line) 
are also presented. 
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Figure 6. Evolutions of V-K,V-J,V-I and B-V colours with the KOI IMF for Models with- (black, line) and without- TPAGB 
(black, line-f solid circles). Also shown are the results with the KOI IMF from the SB99 code, which employs the Padova-AGB (Model 
SB99P-AGB, green, line) and Padova-STD (Model SB99P-STD, green, lino-f-solid circles) evolutionary libraries, respectively. 
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Figure 10. Solar-mctallicity stellar evolutionary tracks (grey line) and isochrones (red line) in the SB99 models. The tracks and isochrones 
(i.e. SB99P-STD and SB99P-AGB) based on the Padova-STD and Padova-AGB libraries are presented in the left and right panels, 
respectively. The masses of stars are in the range of 0.15-120 Mq. The ages of isochrones are log(t/yr)=6, 7, 8, 9 and 10. 
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Figure 11. Left panel is for the comparison in the solar-metallicity isochrone among different authors at ages of log(t/yr)=6, 7, 8, 9 and 
10. Black solid, red dashed and green dot-dashed lines are for the BBCFN-1994, Girardi-2000 and Marigo-2008 [without log(i/yr)=6] 
isochrones, respectively. Blue dot-dot-dot-dashed line, blue triangles and grey dotted line are for the SB99P-AGB, SB99P-STD and 
Yunnan-III isochrones. For the sake of clarity, in the right panel we enlarge the low-T(,fi and high-L region (the region surrounded by the 
grey lines in the left panel). In the right panel, the tracks beyond TP-AGB phase have been neglected for the BBCFN-1994, Girardi-2000, 
Marigo-2008 and SB99P-AGB isochrones. 



4.1 Stellar evolutionary tracks and their 
simplification 

4.1.1 Stellar evolutionary tracks 

We use the MESA version 3709 and the set of input physics 
and parameters described in Section [2] to calculate stellar 
evolution models for a set of solar-metallicity stars (i.e. hy- 
drogen abundance X=0.70 and metallicity Z=0.02). This 
set of stellar evolution models, which are presented by grey 
lines in Fig. (TJ forms the stellar evolutionary library of the 
Yunnan-III EPS models. The masses of stars are from 0.1 
to 100 M© and at a logarithmic interval of Alog(M/MQ) = 
0.05. Using these tracks, we can build the isochrones of SPs. 

TP-AGB stars are obtained by the MESA stellar evo- 
lution code in the Yunnan-III EPS models. This is quiet 



different from the methods of using a synthetic code or 
perturbation function, by which TP-AGB evolution is ob- 
tained and used by the other EPS models. Moreover, TP- 
AGB stars give more contribution to the IR flux. Therefore, 
we give some detailed descriptions of the results concerning 
TP-AGB stars and discussions in Appendix [B] 



4.1.2 Simplification 

In order to conveniently interpolate evolutionary tracks 
and build isochrones (see below), we simplify the stellar 
evolutionary tracks on TP-AGB phase, which is shown 
in Fig. E] [gives the TP-AGB tracks for stars with 1.0 
< A//Mq < 7.9 on the Hertzsprung- Russell (HR) diagram], 
for intermediate- and low-mass stars: we firstly obtain the 
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Figure 3. Evolutions of luminosity L (black line, left y-axis) 
and temperature T^g [green line, right y-axis, /{T^ft) = (logTcff- 
3.2)*8.] for a M = 1.99 Mg star on the TP-AGB phase. Also 
shown arc the average L' and T^^ within a pulse (sixangles) and 
the interpolated points (uparrows, only for luminosity). 
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Figure 7. The differences in the V - ii", V - J, V - / and B-V 
colours (from top to bottom) between Models with-TPAGB and 
without- TPAGB for the KOI (rectangles) IMF. 
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Figure 5. Evolution of TP-AGB stars on the Ig (Toff). vs. Ig (g) 
plane. Also shown are the grid covered by the BaSeL stellar spec- 
tra library at solar mctallicity. 




— with TPAGB /KOI 

— without TPAGB/K01 
logCi/y) = 7, 8, 9, 10 



10'* 2x10^ 
wavelength(A) 



3x10^ 



Figure 8. The fluxes for the with-TPAGB-KOl (black line) and 
without- TPAGB-KOl models at ages of log(t/yr) = 7, 8, 9 and 
10. 



average time t' , luminosity L' and effective temperature T^'ff 

r xdt 

within each pulse by using X = {X = L, TcS or t), 

then obtain the simplified tracks by interpolating between 
these average points. In Fig. [3l we illustrate this process 
as an example of a M = 1.99 Mq star: the blue and pur- 
ple sixangles denote the average L' and T^'ff at epoch t' , 
the purple uparrows are those simplified values by interpo- 
lation. In Fig. [T] the blue points are those simplified tracks 
on TP-AGB phase for different stars. 

4.2 Interpolation of evolutionary tracks and 
construction of isochrones 

One of the important steps in EPS models is to interpolate 
stellar evolutionary tracks by using those tracks computed 
by the MESA code. This step will directly affect the results 



of EPS models. To check this step, we will compare the in- 
terpolated with model tracks of stars. First, we choose an 
intermediate-mass star (M=5.62M0) and a low-mass star 
(A'/=1.78 Mq) from the evolutionary library, then obtain 
their interpolated tracks by using their own two adjacent 
model tracks. In Fig. 21 we give the comparisons between 
the interpolated and model tracks. From it, we see that the 
interpolated tracks are in good agreement with the model 
ones on all phases (including the loop) for both intermediate- 
and low- mass stars. 

Using these model and interpolated evolutionary tracks, 
we build the isochrones of SPs. The ages of SPs are in the 
range of lMyr-15Gyr and at a logarithmic age interval of 
A log(t/yr) = 0.05. In Fig.[Tl the constructed isochrones are 
also presented at ages of log(t/yr) = 6, 7, 8, 9 and 10. 
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Figure 9. The ratio of flux between the without-TPAGB-KOl 
(^A, without) and with-TPAGB-KOl (i^A,with) models at an age of 
log(i/yr) = 8.6. 



4.3 Colours and ISEDs of SPs 

Combining with BaSeL-2.2 stellar spectra library and the 
KOI IMF, we present the colours and ISEDs of solar- 
metallicity SPs in the range of lMyr-15 Gyr at a logarithmic 
age interval of 0.05 (84 ages). The S55 IMF is the most com- 
monly used IMF, nearly all EPS models present the results 
via the S55 IMF. Therefore, we also present the results for 
models with the S55 IMF. As said in Section 2.2, the BaSeL 
library covers the region spanned by TP- AGE stars on the 
log Teff-log g plane at solar metallcity (see Fig. [SJ , but does 
not comprise TP- AGE stars. That is to say, we use the nor- 
mal M-giant spectra for TP- AGE stars. In Fig. [S] we present 
the evolution of TP- AGE stars computed via the MESA code 
and the grid covered by the EaSeL library at solar metallic- 
ity on the log Teff-log g plane. 

In Fig. |6l we present the evolutions oi V — K, V — J, 
V - I and B - V colours for models with the KOI IMF. 
Eecause that the effect of IMF on the results is small (the 
differences in the colours between the adoption of the S55 
and KOI IMFs are less than one tenth of a magnitude), 
this conclusion also has been made by M98 by using several 
IMFs, in Fig. [6] the results with the S55 IMF are not been 
given. As said above, one main characteristic of the Yunnan- 
III EPS models is that the employed MESA code can evolve 
stellar models through TP- AGE phase for intermediate- and 
low-mass stars. Therefore, in Fig. [6] we also present the re- 
sults when neglecting TP- AGE stars (i.e. neglecting the blue 
points in Fig. [T]). We refer to the results for the SPs with 
and without TP-AGE stars as Models with-TPAGE and 
without-TPAGE, respectively. To distinguish the IMF, the 
name of IMF is the supplement to the name of model. 

From Fig. |S] we see that the inclusion of TP- AGE stars 
can make V — K , V — J and V — R colours of SPs redder 
at log(f/yr) > 7.6, as known from M98 and M05, and the 
longer is the passband-wavelength of one magnitude in a 
colour, the more significant is the discrepancy in the colour 
[for example, d{V — K) > A{V — /)]. The discrepancy in the 
B — V colour of SPs is insignificant at all ages. The maximal 
differences in the V — K, V — J and V — R colours of SPs 
are at log(t/yr)'^8.6 (see Fig. [7]). In Fig. [71 we present the 



differences in the V — K,V — J,V — I and B — V colours 
between Models with-TPAGE and without-TPAGE for the 
KOI IMF. The reason for the above fact is that some stars 
along the isochrone are on TP- AGE phase at log(t/yr) > 7.6 
(see Fig. [TJ corresponding to the onset time of TP-AGE 
phase for a ~8M0 star [see Table iBll and Fig. |Bl]). These 
TP-AGE stars are cooler, and would affect the IR flux but 
almost do not affect the visible flux. 

In Fig. [5] we present the fluxes {Fx) for the with- 
TPAGE-KOl and without-TPAGE-KOl models at ages of 
log(t/yr) = 7, 8, 9 and 10. From it, we see that the in- 
clusion of TP- AGE stars would affect significantly the IR 
flux at ages of log(t/yr) = 8.0 and 9.0. For the sake of 
clarity, in Fig. [9] we present the ratio of flux between 
the without-TPAGE-KOl and with-TPAGE-KOl models 
{Fx 

.without /-Fa, with) at an age of log(t/yr)=8.6. From it, we 
see that Fa, without is ~50 percent of -Fa, with at wavelength 
A ~ 2 ^m (i.e. K-band) for a SP at an age of log(t/yr) ~ 8.6. 

4.4 Comparisons in colours with SB99 models: 

We also obtain the results from the SB99 code, which em- 
ploys the same stellar spectra library (EaSeL spectral li- 
brary) and IMFs (i.e. the S55 with a = -2.35 and KOI 
IMFs) as in this work but different stellar evolution library. 
The SB99 code only uses the solar-metallicity Padova-AGE 
and Padova-STD libraries (i.e. the 44th and 34th tracks, re- 
spectively) in this work. The adopted ingredients of the SB99 
code are listed in Tabled We refer to the SB99 models em- 
ploying the Padova-AGE and Padova-STD libraries as Mod- 
els SE99P-AGE and SE99P-STD, respectively. The prop- 
erties of Padova-STD and Padova-AGE tracks have been 
described in Section |3l The upper and lower mass limits 
of IMFs are 100. and 0.1 Mq. For the SB99-v6.02 code, we 
choose such a set of input parameters and physics for the 
following reasons: discuss the sole effect of stellar evolution 
and do not introduce the effects of stellar spectra library 
and IMF on the results. 

In Fig. |6] we also give the results with the KOI IMF 
from the SB99 code, do not present the results with the S55 
IMF from the SB99 code because that the effect of IMF 
on the colours of Models SE99P-AGE and SE99P-STD is 
small. In the following, the results from the SB99 code and 
Yunnan-III models, unless otherwise indicated, mean those 
with the KOI IMF. From Fig. [g] we see that the V - K, 
V - J and V - I colours of Models with-TPAGE/without- 
TPAGE are very different from those from the SB99 code 
for 6.6 < log(t/yr) < 7.6. This is caused by the large dis- 
crepancy in the stellar evolutionary tracks of massive stars. 
When 6.6 < log(t/yr) ^7.6, the IR fluxes of SPs mainly are 
contributed by massive stars ( ^ 8 M© at solar metallicity in 
this work), which evolutionary results are prone to be af- 
fected by the input parameters. In Fig. 1101 we present the 
stellar evolutionary tracks and isochrones of the SB99 mod- 
els. The isochrones are obtained with the help of the SB99 
code and built on the basis of the Padova-STD and Padova- 
AGE libraries, respectively. From it, we see that the massive 
stars in the Padova-AGE and Padova-STD libraries evolve 
through WR phase, and do not extend to very cool region in 
comparison with those in our stellar evolution library (see 
Fig. 1). 

In general, at log(t/yr) > 7.6, the colour-evolution 
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trends of Models with-TPAGB/without-TPAGB agree to 
Models SB99P-AGB/SB99P-STD, respectively. That is 
to say, the isochrones of Models with-TPAGB/without- 
TPAGB coincide with those of Models SB99P-AGB/SB99P- 
STD, respectively. In Fig. Illl we give the Yunnan-III, 
SB99P-AGB and SB99P-STD isochrones at ages of log(i/yr) 
= 6, 7, 8, 9 and 10. From them, we indeed see that the 
Yunnan-III isochrones agree to the SB99P-AGB/SB99P- 
STD isochrones. To be exactly, at log(f/yr) > 7.6, there exist 
differences in the colours between Models with-TPAGB and 
SB99P-AGB, and those between Models without-TPAGB 
and SB99P-STD. They are as follows. 

(i) The colours (except V — I) of Models with- 
TPAGB/without-TPAGB are smaller than the correspond- 
ing ones of Models SB99P-AGB/SB99P-STD by less than 
0.1 mag at large ages [log(f/yr) ^9.0]. 

(ii) The colour evolution of Model with-TPAGB is 
smoother than that of Model SB99P-AGB in the range 
7.6 < log(t/yr) < 8.5, there is a jump in the V — K, 
V -^J and V - /"colours for Model SB99P-AGB at age 
log(i/yr)~8.0. Within 7.6 < log(t/yr) < 8.0 the colours of 
Model with-TPAGB are larger (~0. 2-0.4 mag), but within 
8.0 < log(t/yr) < 8.5 are smaller (~0. 2-0. 4 mag) than those 
of Model SB99P-AGB. All above phenomena are partly due 
to the difference in the maximal mass of stars experiencing 
TP-AGB phase M,nax,TPAGB. In the Padova-AGB library, 
Mmax.TPAGB is 5 Mq , whilc lu our stellar evolution library 
it is greater than 8 Mq (refer to Table |BT] and Fig. [BTI 
in Appendix IB1|I at solar metallicity. Therefore, the jump 
in the colours of Model SB99P-AGB appears at log(t/yr) 
~ 8.0, this epoch corresponds to the TP-AGB phase of a 
M = Mmax.TPAGB (i-c. 5 Mq ) Star in the Padova-AGB li- 
brary (see the 5l"^ line for a Z—0.02 and 5 Mq star in the 
Padova-AGB library from the SB99 code). In the range of 
7.6<log(t/yr) <8.0, there are no TP-AGB stars along the 
SB99P-AGB isochrones, but there are TP-AGB stars along 
the Yunnan-III isochrones. Therefore, the colours of Model 
with-TPAGB are redder than those of Model SB99P-AGB. 
In the range of 8.0 < log(t/yr) < 8.5, there are TP-AGB 
stars along the SB99P-AGB isochrones, and their TP-AGB 
stars are cooler and redder than those on the Yunnan-III 
isochrones (see Figs. [TOl and [TTI only 8.0 < log(t/yr) < 9.0). 
Therefore, the colours of Model with-TPAGB are bluer than 
those of Model SB99P-AGB. 

(iii) The contribution of TP-AGB stars to the IR flux 
and colours peaks at log(f/yr)=8.6 for the Yunnan-III mod- 
els. This epoch corresponds to the TP-AGB phase for 
a Z=0.02, M ~3.16M0 star (refer to Table ^ in Ap- 
pendix lBTj) . However, the maximal contribution of TP-AGB 
stars is at log(t/yr)=8.0 for the SB99 models with the Padova 
tracks, this corresponds to the epoch of jump in the colours. 



5 COMPARISONS WITH THE OTHER EPS 
MODELS 

In this section we will compare our results with the other 
studies, which have considered TP-AGB stars in their EPS 
models. For these models, their main characteristics have 
been described in Section O 

Here, we do not use all the models, only use the CB07, 
M05, VIO and POPSTAR models at solar metallicity. The 



comparisons in the colours with the SB99 models have been 
presented in the Section 14.41 (i) For the CB07, we use the 
set of models by using the Padova-1994 stellar evolution- 
ary tracks, BaSeL-3.1 stellar spectra library, the S55 and 
Cha03 IMFs (which we call CB07-S55 and CB07-Cha03). 
(ii) For the POPSTAR, we use the set of results by using 
the S55 with the lower and upper mass limits of 0.15 and 
100 M0, KOI and Cha03 IMFs (which we call POPSTAR- 
S55, POPSTAR-KOl and POPSTAR-Cha03). (iu) For the 
VIO, we choose the set of models by using the KOI and 
KOl-cor IMFs (which we call VlO-KOl and VlO-KOl-cor). 
(iv) For the M05, we use the set of models by using the S55 
and KOI IMFs (which we call M05-S55 and M05-K01). The 
corresponding ingredients for the used EPS models are given 
in TableO It is worth mentioned that at solar metallicity the 
M05 and VIO models indeed use the 'Cassisi-|-Geneva' tracks 
and GOO isochrones, respectively (in comparison with Ta- 
ble 1). Moreover, the SB99P-AGB/SB99P-STD isochrones 
corresponds to the Marigo-2008/Girardi-2002 isochrones, re- 
spectively, because the same evolutionary tracks but differ- 
ent construction method are used. 

We choose the above EPS models for the following rea- 
sons: keep the common ingredients as more as possible. From 
Tabled we see that these EPS models have one or several 
ingredients in common. This guarantees that the sole effect 
of ingredient on the results can be derived. M05 has men- 
tioned that EPS models keep the uncertainties inherent in 
the stellar evolutionary tracks and in the spectral transfor- 
mations, the isolation of the effect on the final results is 
important and quantifying the uncertainties is very impor- 
tant since the cosmological inferences that are derived on 
the basis of galaxy ages and metallicities ultimately rely on 
the SP models. 



5.1 Comparison in the isochrone 

For the sake of the paper's construct, we first present the 
comparison in the isochrone among these EPS models be- 
fore analysing the differences in the results. The POPSTAR 
models used the Padova-1994 tracks (together with the re- 
vision of the BGS98 isochrone, see Table [l}, from which 
the BBCFN-1994 isochrones are derived (see Section O. 
At solar metallicity, the VIO models indeed used the GOO 
isochrones (see Table [1} , which are a part of the Girardi- 
2002 isochrones (see Section[3|. The Padova-1994 tracks and 
the prescriptions of MG07 in the used CB07 models (see the 
first paragraph of Section [5} correspond to the Marigo-2008 
isochrone (see Table [T] and Section [3] for details). Therefore, 
in Fig. [ni we also present the BBCFN-1994, Girardi-2002 
and Marigo-2008 isochrones at ages of log(f/yr)= 7, 8, 9 and 
10 and Z=0.02. These three sets of isochrones are from the 
Padova's website. For the BBCFN-1994 isochrones, we use 
Mbol=4.83. 



5.2 Comparisons in the colours 

In Fig. [HI we present the V - K, V - J, V - I and 
B - V colours for the with-TPAGB-KOl, SB99P-AGB- 
KOl, CB07-S55, CB07-Cha03, M05-S55, M05-K01, VlO- 
KOl, VlO-KOl-cor, POPSTAR-S55, POPSTAR-KOl and 
POPSTAR-Cha03 models, (i) From the lower right panel 
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Figure 12. Comparisons in tlie V — K, V — J , V — I and B — V colours among the with-TPAGB-KOl (red open rectangles), SB99P- 
AGB-KOl (green open triangles), CB07 (blue, solid line for CB07-S55, solid line+circles for CB07-Clia03), VIO (cyan, dashed line for 
VlO-KOl-cor, dashed line+circles for VlO-KOl), M05 (magenta, dot-dashed line for M05-S55, dot-dashed line+circles for M05-K01) and 
POPSTAR (black, dotted line for POPSTAR-S55, dotted line+circles for POPSTAR-KOl, dotted line+crosses for POPSTAR-Cha03) 
models. 



Table 2. Main characteristics of the used EPS {SB99, POPSTAR, VIO, M05 and CB07) models. The 2nd and 3st columns are for 
mctallicity and age range. The 4th and 5th columns are library of stellar evolution and the name of the corresponding isochrone. The 
6th and 7th columns arc spectral library and IMFs. 



Model 


Z 


Age 


stellar evolution 


corresponding 


stellar spectra 


IMF 






log(Vyr) 


library 


isochrone 


library 




Yunnan-III 


0.02 


6.0-10.18 


Yunnan-III 


Yunnan-III 


BaSeL 


S55/K01 


SB99 


0.02 


6.0-10.18 


Padova-AGB [MG07] 


SB99P-AGB 


BaSeL 


S55/K01 








Padova-STD [GOO] 


SB99P-STD 






POPSTAR 


0.02 


5.0-10.18 


Padova-1994+BGS98 


BBCFN-1994 


BaScL+othcr 


S55/K01/Cha03 


VIO 


0.02 


7.8-10.25 


GOO 


Girardi-2002 


MILES 


/KOl/KOl-cor 


CB07 


0.02 


5.0-10.30 


Padova-1994+MG07 


Marigo-2008 


BaScL-3.1 


S55/ /Cha03 


M05 


0.02 


3.5-10.18 


Cassisi + Geneva 




BaSeL-3.1+LM02 


S55/K01 



of Fig. 1121 we see that the B — V colour agrees with each 
other (except VlO-KOl-cor) when Iog(t/yr) > 6.6. (ii) From 
the lower left panel of Fig. 1121 we see that the V — I 
colour is in agreement for all models (except VlO-KOl- 
cor) when log(t/yr) > 6.6, but is not so good as that of 
B — V colour, the discrepancy in the V — I colour among 
all models is within ~0.5mag. When 8.4 < log(f/yr) < 9.0 
and log(i/yr) > 9.0, the V - J colour of MOslmd POPSTAR 
models is greater than the other models, respectively, (iii) 



From the top-right panel of Fig. 1121 we see that the V — J 
colour coincides for all models (except VlO-KOl-cor) when 
log(t/yr) >6.6. But, at smaU [6.6 < log(t/yr) < 7.2] and in- 
termediate [8.0 < log(t/yr) < 9.6] ages, the V — J colour 
of the with-TPAGB-KOl models (< 0.6 mag) and those of 
the CB07 and M05 models ( < 0.6 mag) are significantly 
greater than the value of the other models, respectively, 
(iv) From the top-left panel of Fig. 1121 we see that the 
V ~ K colour agrees with each other (except VlO-KOl-cor). 
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However, at small ages [6.6 < log(t/yr) < 7.2], the V — K 
colour of with-TPAGB-KOl and M05 models is greater 
(< 1.0 mag) than the other models. At intermediate ages 
[8jO<log(t/yr) <9.6], the V - K colour of CB07 and M05 
models is significantly greater (< 1.0 mag) than the other 
models. The following, we will give the comparisons in the 
colours between ours and other models and analyse the rea- 
sons. 



5.2.1 Comparison between imth-TPAGB and POPSTAR 
models: 

The colours of POPSTAR models agree with those of with- 
TPAGB-KOl models except 1/ - A", 1/ - J and V - / in the 
range of 6.6 < log(t/yr) < 7.6. The colour of POPSTAR 
models is redder than that of with-TPAGB-KOl models at 
large ages [~0.2mag, log(t/yr) ^9.0]. Moreover, there is no 
jump in the evolution of V — K colour (exists for SB99P- 
AGB models) for POPSTAR models. The effect of IMF (S55 
& KOI & Cha03) on the colours [except V — I colour at 
7.0 < log(t/yr) < 8.0] is insignificant, the differences in the 
colours between the POPSTAR-S55 and POPSTAR-Cha03 
models are small. 

POPSTAR-KOl models use the BaSeL-|-other Ubrary 
and KOI IMF. So the difference in the V — I colour be- 
tween the POPSTAR-KOl and with-TPAGB-KOl models at 
large ages is caused by the differences in the stellar evo- 
lutionary tracks, spectral library and some techniques in 
the construction of EPS models. At large ages, SPs com- 
pose of low-mass stars. In these two sets of stellar evolu- 
tionary tracks used by the POPSTAR and with-TPAGB- 
KOl models, the method obtaining TP-AGB tracks (simply 
synthetic code or complete evolution code), input physics 
and parameters (such as opacity, EOS and overshooting) 
are different. The stellar evolutionary tracks used by the 
POPSTAR and with-TPAGB-KOl models correspond to 
the BBCFN-1994 and Yunnan-III isochrones (see Table [2]). 
From Fig. [Til we see that the first GB along the BBCFN- 
1994 isochrones has lower temperature than that along the 
Yunnan-III isochrones at log(i;/yr) = 10., this leads to red- 
der V — I colour. Because the first GB does not extend to 
very low-temperature region, therefore, it does not affect the 
V — K and V — J colours. 



5.2.2 Comparison between with-TPAGB and VI models: 

The beginning time of the VIO models is log(t/yr)~7.8. 
First, from Fig. \U[ we see that the effect of IMF (KOI & 
KOl-cor) on the colours is significant. The results of the 
VlO-KOl-cor models are quiet different from those of the 
VlO-KOl and the other EPS models. The results of the VIO- 
KOl models agree with those of the with-TPAGB-KOl mod- 
els. However, at intermediate ages [8.0 < log(t/yr) ^9.0], the 

V — K and V — J colours of the VlO-KOl models are redder 
(~0.4 and 0.3 mag); and at large ages [log(t/yr) > 9.0], the 

V — I colour is bluer (~0.1 mag) than the corresponding one 
of the with-TPAGB-KOl models. 

The VlO-KOl models use the GOO isochrones, which is 
a part of the Girardi-2002 isochrones. Moreover, the VlO- 
KOl models use MILES stellar spectra library (MILES com- 
prises some TP-AGB stars although at non-solar metallic- 



ity), but the with-TPAGB-KOl models use the BaSeL li- 
brary. Therefore, the differences in the colours between the 
VlO-KOl and with-TPAGB-KOl models are mainly caused 
by the adoptions of different spectral and stellar evolution 
libraries and different constructional techniques. Compar- 
ing the isochrones in Fig. 1111 we see that the TP-AGB stars 
along the Girardi-2002 isochrones do not extend into the 
very cool region than those along the Yunnan-III isochrones 
when 8.0 < log(t/yr) ^9.0, and the first GB has lower tem- 
perature at large ages, the V — K and V — J colours at in- 
termediate age should be bluer and the V — I colour at large 
ages should be redder. But the fact is the contrary. There- 
fore, the difference in the stellar spectral library (MILES) 
should be the main factor that causing the differences be- 
tween the VlO-KOl and with-TPAGB-KOl models. 

Recall the difference in the V — K and V — J colours 
between the SB99P-AGB and SB99P-STD models in Fig. [6] 
and the fact that the SB99P-STD isochrones correspond to 
the Girardi-2002 isochrones (the VIO and SB99P-STD mod- 
els use the same stellar evolutionary tracks at intermediate 
and large ages), we conclude that the V — K and V — J 
colours will be redder by ~ 0.5 mag if the VlO-KOl models 
use the Marigo-2008 isochrones, also true for the Yunnan-III 
isochrones. This is because that the Yunnan-III isochrones 
are similar to the SB99P-AGB isochrones (conclusion made 
in Section [4.4|) and the SB99P-AGB isochrones correspond 
to the Marigo-2008 isochrones. 



5.2.3 Comparison between with-TPAGB and M05 models: 

First, from Fig. [12] we see that the effect of IMF (S55 & 
KOI) on the colours is insignificant except V — K and V ~ ,J 
colours in the age of 7.6 < log(t/yr) < 8.4. The V-K,V-,J, 
V — I and B — V colours are redder than the corresponding 
one of the with-TPAGB-KOl models by ~lmag, ~0.6mag, 
^0.3 mag and ~0.1mag in the age of 8.2 < log(t/yr) ^9.0. 

The M05 models use the BaSeL library, LM02 spec- 
tra for TP-AGB stars and the 'fuel consumption' theorem, 
therefore, the differences in the colours between the M05- 
KOl and with-TPAGB-KOl models are partly introduced by 
the adoption of 'fuel consumption' theorem and the spectra 
of TP-AGB stars. 



5.2.4 Comparison between with-TPAGB and CB07 
models: 

The colours of the CB07 models agree with those of with- 
TPAGB-KOl models at log(t/yr) > 6.6 except V - K &nd 
V — J colours at intermediate ages [redder by --^1.0 mag, 
8.0 < log(f/yr) ^9.6] and V — I colour at large ages [bluer 
by ~0.2mag, log(t/yr) > 9.0, agrees with VIO]. The effect of 
IMF (S55 & Cha03) on the colours is insignificant except at 
large ages [log(t/yr) > 10.0, a hundredth of a mag] for the 
CB07 models. 

The effect of IMF on the colours is small for both CB07 
and with-TPAGB models, and two models give the results 
with the S55 IMF. CB07 models use BaSeL spectral library 
and the Padova-1994 tracks together with the description of 
MG07 (corresponds to the Marigo-2008 isochrones). There- 
fore, The differences, in the V — K and V — J colours at in- 
termediate ages and V — I colour at large ages, between the 
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CB07-S55 and with-TPAGB-S55 models are caused by the 
difference in the stellar evolution tracks and the techniques 
during the construction of EPS models. From Fig. 1111 we see 
that TP-AGB stars on the Marigo-2008 isochrones extend 
into cooler region than those on the Yunnan-III isochrones 
at intermediate ages, therefore, the redder V — K and V — J 
colour. Moreover, from Fig. 1111 we indeed see that the tech- 
nique of building isochrones would be a reason, the Marigo- 
2008 and SB99P-AGB isochrones use the same evolutionary 
tracks, but the isochrones are different. 

5.^2.5 Comparison among all models: 

Within 6.6<log(f/yr)<8.0, the SB99P-AGB, CB07 and 
POPSTAR models (KOl-IMF better) are in agreement with 
each other, the VIO models do not cover this age range. 
The main reason for the the agreement is that they use the 
same stellar evolutionary tracks for massive stars (i.e. the 
Padova-1994, see Tabled) and stellar spectra library. The 

V - K (except for M05), V - J and V - I colours of the 
with-TPAGB-KOl models are redder than all of the other 
models. The reason for the differences has been discussed in 
Section 1331 

At large ages log(f/yr) > 9.0, all colours of all models 
(except VlO-KOl-cor) are in agreement. The CB07 mod- 
els have reddest V - K and V - J colours; POPSTAR has 
the reddest V - I colour (> with-TPAGB/SB99P-AGB > 
CB07/V10 > M05 in turn). 

At intermediate ages [8.0 < log(t/yr) < 9.0], the V - K 
and V — J colours of CB07 and M05 models and the V — I 
colour of M05 model are redder than all the other models. 
Moreover, the maximal disagreement in the colours among 
different literatures appears in this range, especially in the 

V — J and V — K colours. 

5.3 Comparison in the contribution of TP-AGB 

stars 

In Section 21 we conclude that the maximal contribution 
of TP-AGB stars to the /T-band is about 50 per cent. 
This value i s some what different from the other authors: 
iFrogel et al] (|l990l ) concluded that the TP-AGB contribu- 
tion to the total luminosity of SP reaches a maximum of 
about 40 percent at ages from 1 to 3 Gyr, falling to less 
thar i 10 perc e nt at 10 Gyr (from M08); 

iBruzuall (|201lh mentioned that it can reach to 60 per- 
cent in the CBll models, close to a factor of two more than 
in the BC03 model; M05 found that 40 percent of the bolo- 
metric contribution and 80 percent in K passband. Part of 
reasons is that we used the BaSeL library, the spectra of 
TP-AGB stars are represented by normal M-type giants. 
Another reason is the less redder TP-AGB stars and shorter 
during-time are produced when comparing with some EPS 
models. 



6 SUMMARY AND CONCLUSIONS 

Using the MESA stellar evolution code, we compute the con- 
sistent stellar models for low-, intermediate- mass and mas- 
sive stars in the range of 0.01-100 M©. By combining with 
BaSeL stellar spectra library and the KOI and S55 IMFs, 



we build the Yunnan-III EPS models, and present colours 
and ISEDs of solar-metallicity SPs in the range of IMyr- 
15 Gyr. By comparing the results between Models with- 
TPAGB and without-TPAGB, we confirm that the inclusion 
of TP-AGB stars can raise V — K, V — R colours and IR 
fiux for SPs at log(t/yr) > 7.6, the maximal value is at an 
age of log(t/yr)~8.6 (about 0.5-0.2 mag for colours, about 
2 times for A'— band fiux). Moreover, by comparisons, we 
find that the colour-evolution trends of Model with-TPAGB 
are similar to those obtained by the SB99 code employ- 
ing the Padova-AGB library. At last, by comparing the re- 
sults among ours, CB07, M05, VIO and POPSTAR, we find 
that the difference in the B — V colour is small, but the 
difference in the V — K colour is large (~1 mag when 8. 
< log(t/yr) < 9.) The above conclusion is independent of the 
adoption of^IMF (KOI & S55 IMFs). 

The stellar evolutionary tracks, isochrones, 
colours and ISEDs of SPs can be obtained on re- 
quest from the first author or from our website 
(I http:/ /wwwl.ynao.ac.cn/^zhangfh/[ ). You can build 
your EPS models by combining the isochrones with the 
stellar spectra library and IMF you wanted. Now the format 
of stellar evolutionary tracks we provided is the same as 
that in the SB99 code, you can input easily them in it 
to get the SPs' results. Moreover, the colours involving 
other passbands or on other systems (for example, HST 
F439M^ — F555W colour on AB system) can also be 
obtained on request. In this paper we have only presented 
the results of solar-metallicity SPs - more detailed studies 
will be given later. 
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Figure Bl. The onset time [log(txpAGB,b)! left y-axis] and dur- 
ing [txpagBi right y-axis] of TP-AGB phase as a function of 
stellar mass. The results, obtain by using the standard or the 
C13 set of input parameters and physics in the MESA code, are 
represented by Line+circles and stars, respectively. 
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Figure B2. The relation between T^g and mass loss M for TP- 
AGB phase of stars with 1.0<M/Mq <7.9. The symbols have 
the same meaning as in Fig. [2] These results are obtained by 
using the standard set of input parameters and physics in the 
MESA code. 

APPENDIX A: THE MAGNITUDES AND 
COLOURS FOR MODEL WITH-TPAGB. 

In this section, we present the integrated magnitudes and 
colours of solar-metallicity SPs when considering TP-AGB 
evolution (i.e. Model with- TPAGB). Tables W\\ and IX2l are 
the results when using the KOI and S55 IMFs, respectively. 



APPENDIX B: EVOLUTIONARY RESULTS OF 
TP-AGB STARS 

In this section, we will present the evolutionary results of 
TP-AGB stars. In Appendix IBll we present the onset time, 
lifetime and inter-pulse time of TP-AGB phase, the maximal 
mass of star experiencing TP-AGB phase, and some com- 
parisons with those of MG07. In Appendix IBll all results are 
obtained by using the set of input physics and parameters, 
which have been described in Section [2] and called the stan- 



Table Bl. The onset time tTPAGB.b (the 2nd column), duration 
TTPAGB (the 3st column) and the pulsation time nxpAGB (the 
4th column) of TP-AGB phase for solar-metallicity stars with 
different masses (the 1st column). Top panel is for the standard 
set of parameters and physics in the MESA code. Bottom panel 
is for the C13 set of parameters and physics (which is described 
in Appendix IB2I I. 



Mass 


*TPAGB,b 


TTPAGB 


'T'TPAGB 


(Mo) 


(Myr) 


(Myr) 




standard set 


1.00 
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0.2549 
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0.8362 
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14 


2.00 
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1.4787 
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14 


C13 set 


1.12 


7854.9492 
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1.41 
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0.8694 


11 
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18 


4.47 
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5.62 
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0.1219 


6 


7.08 


49.7702 
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14 



dard set, in the MESA code. In Appendix lB2l we present the 
results by using another set of input parameters and physics 
used in the MESA code, and discuss the effects on the onset 
time and lifetime of TP-AGB phase. 

Bl Some results by using the standard set of 
physics and parameters 

Bl.l Onset time and duration of TP-AGB phase: 

First, in Table lBll we present the onset time (fxpAGB.b, accu- 
rately speaking, it is the average time during the first ther- 
mal pulse t'l), the pulsation time (titpagb) and the duration 
or lifetime (ttpagb, the difference between the average time 
during the last thermal pulse iuTPACB ^"-^ '^^ TP-AGB 
phase as a function of stellar mass. Then, in Fig. IBll we 
give log(tTPAGB,b) and ttpagb as a function of stellar mass, 
which is in the range of 1.0-7.9Mq. 

From them, we see that the logarithmic onset-time of 
TP-AGB phase [log(tTPAGB,b] hnearly decreases with the 
logarithmic initial-mass log(i\f). The lifetime of TP-AGB 
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Table Al. The bolometric magnitude Mbol, V-magnitude My, U-B, B-V,V-R, V-I,V-J,V-H and V - K Colours of 
solar-metallicity SPs when considering TP-AGB stars (Model with-TPAGB) and using the KOI IMF. 
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phase reaches the maximum (~ 2Myr) at M—2.24Mq. This 
differs from that of MG07: the majcimal during of ~ 4 Myr at 
M=2Mq for M-type TP-AGB stars (see the left-top panel 
of their Fig. 22). That is to say, the maximal lifetime of TP- 
AGB phase is haff of the MG07 result. 



B1.2 the maximal mass of stars experiencing TP-AGB 
phase: 

Here need to be mentioned is the maximal mass of stars ex- 
periencing TP-AGB phase (i\fmax,TPAGB). MG07 have ob- 



tained that the upper mass limit for the development of the 
AGB phase is located close to 5Mq for all values of Z as a 
consequence of overshooting. This value of 5Mq is smaller 
than th at of 8 Mq found i n e arlier models (see the review 
paper of lvan Winckeill2003l and lHerwidbOOSl ). Moreover, the 
classical results found an important TP-AGB contribution 
at 5 M0 in overshooting tracks (Padova) and at 3 Mq in 
tracks without overshooting (M98). Therefore, 5Mq is the 
highest initial mass to be considered in the MG07 work, but 
they also mentioned that the value of Mmax.TPAGB is slightly 
uncertain because the initial phase of carbon burning has not 
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Similar to Tabic lAll but for the results when using the S55 IMF. 
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M >8M0, the timescale on TP-AGB phase is very short 
and the contribution is smaller. 



B1.3 the minimal temperature for superwind: 

In Fig. IB2I we present the TP-AGB stars on the Toff and 
mass-loss M plane, from it we see that the minimal temper- 
ature for superwind is log(Tcfl /K)~3.5. 



been followed in detail. In the FSPS models (described in 
Section [S]), it is cited from M08 that Mniax,TPAGB is 5M0. 

In our calculations, Mmax.TPAGB is greater than SMp 
at solar metallicity. In the Fig. 22 of iPaxton et al.l (|201ll ). 
they presented a grid of MESA star evolutionary tracks 
with mass ranging from 2 to 10 Mq with Z—0.02, and in 
their paper they said that the 8 and 10 M0 models start 
to ignite carbon burning off center, this implies that these 
stars are experiencing envelope burning which keeps the TP- 
AGB contribution low, where 2-7 Mq models do not and will 
presumably go on form C/0 core white dwarfs. Therefore, 
the Mmax.TPAGB is greater than 8M0, and for stars with 
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B2 Results by using another set of physics and 
parameters 

In this work, we only calculate the stellar evolution models 
and present the SP's results at solar metallicity. This can 
not prove the reliability of the standard set of physics and 
parameters, in the MESA code, on the TP-AGB evolution. 
Therefore, we change some input parameters and physics in 
the st andard set by mimicing the set used bv lPaxton et al.l 
(|201lf ) to explain the formation of a C13 pocket (M =2 
and ^=0.01). This set of physics and parameters is called 
the C13 set. 

In the C13 set, the changed parameters and physics 
are as follows, (i) Overshooting: at both boundaries for H- 
burning and non-burning convective zones, / (see Section 
2.1) is changed from 0.01 to 0.014. Moreover, overshoot- 
ing is allowed for He-burning and Z-burning zones, and at 
both boundaries /=0.014. During the TDU /=0.126 at the 
bottom of the convective envelope and /=0.008 at the bot- 
tom of the H-shell flash convection zone, (ii) Opacity: CO 
enhanced opac ity is used, (iii) Mass loss: the coefficient of 
iBloeckerl (|l995l ) mass- loss-law 77B is changed from 0.1 to 0.05. 
(iv) Atmosphere: the usage of simple photosphere model is 
changed to the usage of model atmosphere tables for pho- 
tosphere. Similarly, we use the C13 set of input parameters 
and physics to obtain the luminosity of 1 Mq star at an age 
of 4.57xl0Vr, log(L/L0)=O.O12 for the GS98 opacity while 
log(L/L0)=O.O12 for the GN93 opacity. These luminosities 
are within one hundredth of a magnitude. 

Using the C13 set of physics and parameter in the MESA 
code, we obtain the onset time fTPAGB.b, lifetime ttpagb 
and inter-pulse time titpagb of TP-AGB phase. These re- 
sults are given in the bottom panel of Table [BT] Moreover, 
in Fig. IBll we also plot log(tTPAGB,b) and ttpagb by using 
the C13 set in the MESA code. By comparisons with the 
results of the standard set, we found that the variation of 
iTPAGB.b is little, the ttpagb of low-mass stars becomes to 
be greater (by ~0.2dex) and ttpagb of intermediate-mass 
stars smaller (no more than 0.05 dex) than the correspond- 
ing one, and the maximal lifetime, obtained by using the 
C13 set in the mesa code, can not reach to the value of 
MG07 (the maximum of ~ 4Myr). The input physics and 
parameters in the work of MG07 are very differ from those in 
the MESA code. In the following we describe the parameters 
used by MG07. 

• O pacity: they u se variable molecular opacity by means 
of the iMarigol (|2002l , hereafter M02) method, i.e. opacities 
are consistent with the changing photospheric chemical com- 
position of AGB envelope. The M02 method is different from 
the commonly used one, i.e. interpolations are preformed as 
a function of temperature, density and hydrogen content for 
a fixed initial metallicity (MG07). 

• Mass loss: MG07 use formalisms for the mass-loss 
rates derived from pulsating dust-driven wind models of C- 
and O-rich AGB stars, and the value of mass-loss depends 
on whether a stellar model is oxygen-rich or carbon-rich. 

• Overshooting: MG07 use two parameters: Ac and Ae. 
Ac describes the extent across the border of the convec- 
tive zone and Ac is the overshooting at the lower bound- 
ary of convective envelope. Both are expressed in units 
of pressure scale height. When M<1M0, M>1.5M0 and 
1<M/M0<1.5, Ac is 0., 0.5 and M-1.0, respectively. In the 



stage of core helium burning Ac = 0.5. When M<O.6M0, 
O.6<M/M0<2.O, 2.O<M/M0<2.5 and M >2.5M0, Ac is 
0., 0.25, 0.25-0.7 and 0.7. 

• Temperature: the T^s of TP-AGB models is derived 
with the aid of complete integration of static envelope 
models, extending from the photosphere down to the core 
(MG07). 

The different choice of all these parameter (opacity, the 
mass loss, overshooting, atmosphere) would alter or affect 
the evolution of TP-AGB stars. For example, the adoption of 
mass-loss law would affect TP-AGB evolution (M05). There- 
fore, there exists disagreement between the MG07 and our 
models. 

Similarly, in the 'fuel consumption' technique, uncer- 
tainties in mass loss, mixing and efficiency of hydrogen burn- 
ing at the bottom of the convective envelope together also 
affect the TP-AGB evolutions (see M05). These disagree- 
ment further affects the SP's results with 'fuel consumption' 
techni que. M98 has made a conclusio n that the adoption 
of the lBloecker fc Schoenbernerl (|l99ll ) mass- loss law would 
make the models luminous and evolution faster than that 
predicted by the standard core mass- luminosity relation for 
7M0 star, resulting drastic decrease of the TP-AGB lifetime 
and smaller contribution of TP-AGB phase to the ISED. 

This paper has been typeset from a T^jiX/ M^gX file prepared 
by the author. 



